Quantum Zeno effect by indirect measurement: The effect of the detector 
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We study the quantum Zeno effect in the case of indirect measurement, where the detector does 
not interact directly with the unstable system. Expanding on the model of Koshino and Shimizu 
[Phys. Rev. Lett., 92, 030401, (2004)] we consider a realistic Hamiltonian for the detector with 
a finite bandwidth. We also take explicitly into account the position, the dimensions and the 
uncertainty in the measurement of the detector. Our results show that the quantum Zeno effect is 
not expected to occur, except for the unphysical case where the detector and the unstable system 
overlap. 



I. INTRODUCTION 

In a recent paper [l| , Koshino and Shimizu (KS) con- 
sidered the quantum Zeno effect (QZE) HHIlll, for 
an exactly exponentially decaying system. They conclude 
that the possibility for observing the QZE exists even in 
this case, where the initial deviation from exponential be- 
haviour, thought to be of vital importance for the QZE, 
is absent. 

As an example, they considered a two-level atom 
(TLA) decaying to its ground state by emitting a pho- 
ton counted by a detector. Through a continuous indi- 
rect measurement of the emitted photon, they obtain the 
QZE even in the extreme case where the "jump time" is 
zero, which leads them to the conclusion that the QZE 
is easier than expected so far to occur. 

Since this contrasts with conventional wisdom, we un- 
dertook a careful reexamination of the problem. We find 
that it is essential to reformulate the Hamiltonian so as 
to account for the influence of the finite extent of the 
detector, including its distance from the TLA. Our cal- 
culations, based on a discretization technique and the 
numerical solution of the resulting system of differential 
equations, show that the QZE does not occur, except for 
the unphysical situation, inherent in the model of Ref. jj] , 
in which the TLA and the detector overlap; i.e. the de- 
tector contains the TLA. 



II. HAMILTONIAN CONSTRUCTION 

The system we consider follows as close as possible the 
lines of Ref. [l| (the same system and formalism has 
been employed by KS earlier in Ref. [g)- The unstable 
system, a two-level atom (TLA) with \g) the ground and 
|e) the excited state, is initially in |e) and decays to its 
ground state by emitting a photon. The emitted pho- 
ton is subsequently detected and the "observer" becomes 



aware of the TLA decay. The total quantum system we 
consider includes, besides the TLA and the electromag- 
netic field, a part of the measuring apparatus, which is 
treated quantum mechanically 

The system Hamiltonian (?i = c = 1) in the form em- 
ployed by KS is: 



7i = Tio + ^1 + ^2 
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U2 = // dkduj (V^^kCkt^ + H.c.) + Locl.^C]^uj (4) 
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where Tio is the part representing the free evolution of the 
TLA, Til the atom-photon interaction and the free evo- 
lution of the photon, with h\^ the annihilation operator 
for the photon of k wave vector. The combined system 
TLq + Hi is coupled to a (macroscopic) detector, a part 
of which is modeled by Ti.2 which represents quantum 
mechanically the measuring procedure, i.e. the detection 
of the emitted photon, ^k and y^ are the atom-photon 
and the photon-detector couplings, respectively. All pho- 
ton modes are coupled with the continuum of the bosonic 
elementary excitations in the detector, with annihilation 
operator cu_uj. The usual commutation relations for the 
6, c operators hold. 

We wish to elaborate on two issues on this form of the 
Hamiltonian. First, in TC2, the detection process is ac- 
complished by transferring a quantum of a photon mode 
to the detector modes through the term (feJ^Ckw) which 
conserves k. This means that there is no uncertainty in 
the detection process of k which is a rather unphysical 
feature. Consider a detector capable of detecting (prac- 
tically) all photons. In the case that the electromagnetic 
field decays inside the detector as e~'^'^ (see Fig. ^|, the 
momentum of the detected photon is determined within 
7 (from the uncertainty relations AxAk ^ 1). Thus, 
there is an inherent uncertainty on the outcome of the 
measurement, a photon with k wavenumber can be de- 
tected as k' inside the bandwidth 7. We take this into 
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FIG. 1: A schematic plot of the TLA and the detector 
(shaded region) as we consider it. The electromagnetic field 
emitted by the TLA enters the detector with a penetration 
depth of 1/7. The atom is displaced by xd and the detector 
spans from to xq. 



account by introducing C{k, k') in 7^2 which becomes; 



n. 



dk'dkdtj 



{C{k,k')bl 



Ck'. 



H.c. 



dkcLo 



and C depends on the details of the electromagnetic field 
attenuation inside the detector. 

The derivation of C can be accomplished via two differ- 
ent pathways. In the first, we consider the macroscopic 
characteristics of the decay of the electromagnetic field 
inside the detector and obtain C phenomenologically. In 
the second, we need to specify the details of the detector 
and we can derive C through this more rigorous approach. 
Both ways lead to the same result for C (physically), i.e. 
a smooth function with finite width, which is actually the 
only important factor in our model. We briefly describe 
both. 

In the phenomenological approach we can assume that 
the electromagnetic field attenuation inside the detector 
depends on two factors, the coupling strength (y/rfk) and 
the density p{x) of the bosonic excitations of the detec- 
tor. The latter is introduced to account for a smooth 
transition at the surface to the bulk density and/or for 
other space dependent particulars of the detector. The 
local attenuation rate of the mod e of the electromagnetic 
field is proportional to \/r]kp{x) and the electromagnetic 
field mode (in one dimension) becomes: 



^ — ikx 



Vk{x) = 



, X > Xq 



_;^g-»fcXg-/;« .JnkP(x)dx ^x<Xo 



where A/" is the normalization factor for the photon eigen- 
mode, and xq is the point where the detector begins (see 
FigP). 

The coupling of the electromagnetic field modes (j)k 
whith the detector modes has to be such that their de- 
cay inside the detector is of the form oi Vk- A similar 
approach is followed in Ref.jg] in the context of absorb- 
ing boundaries in spectral methods, where it is shown 
that the couplings C{k,k'), there the coefficients of the 
absorbing boundary linear transform, are the projection 
coefficients of Vk on (pk, where Vk = Vk — (f>k, i-S. the 



part of the Vk mode transferred to the detector. In the 
following numerical analyses, we assume a Gaussian at- 
tenuation inside the detector, which leads in a Gaussian 

C{k,k')^C{k-k'). 

In the case one wishes to take into acount all the details 
of the detector in a more fundamental level the Hamilto- 
nian of the detector and the resulting eigenmodes have 
to be specified. Then the coupling whith the 0^ modes 
is {4>k'\'D\^kuj) where ^kuj are the eigenmodes of the de- 
tector and V the coupling operator of the detector with 
the electromagnetic field. We have to bear in mind that 
the eigenmodes of the detector are spatially localised, i.e. 
inside the detector. Also, since we wish to represent a de- 
tector and not a mirror the eigenmodes have to attenuate 
smoothly at the surface of the detector. Clearly the exact 
calculation has to proceed by an explicit formulation of 
the detector Hamiltonian and determination of its eigen- 
modes. We do not intend to proceed in this direction 
^ina^pur scope in this paper is to demonstrate the qual- 
fetivc effects of the detector width and position of the 
obervation of QZE. The basic result of such an analysis 
can be deduced by considerind a simple form for the $few 
in conformation with the two restrictions we mentioned: 
space localisation and smooth variations, for example a 
plane wave with a Gaussian envelope. In this simple case 
it is evident that the C{k,k') could practically be thought 
of as a Gaussian. 

The 7^2 in T] is a subclass of this generalized version 
with C{k,k') being a delta function. In retrospect, this 
means that 7-^-0, which implies that the physical di- 
mensions of the detector tend to infinity. The latter is a 
direct artificial influence on the dynamics of the decaying 
two-level atom, since it implies that the (inflnite) detec- 
tor and the TLA overlap. We return to this issue latter 
on. 

The second issue we wish to take into account is the 
relative position of the detector and the TLA. This is 
straightforward, and is accomplished by including the 
correct displacement phase factor in the Hamiltonian. 
This phase has the simple form (pk = e*'^^" (see Fig. ^|, 
as employed for example in |3 for the somewhat similar 
case of a TLA coupled through the electromagnetic field 
with another TLA. The way the Hamiltonian is written 
so far, the TLA and the detector overlap and we have to 
displace one of them. It is more convenient to displace 
the atom, since it involves inclusion of the phase factor in 
fewer terms, so the term Hi of the Hamiltonian becomes: 



Hi 



dk 



(Cke''"° |e) (5|6k + H.c.) + fc^J^foJ ■ (6) 



In general the displacement is determined by the prob- 
lem at hand, but in all cases it should be such that the 
atom does not overlap with the detector. In the present 
form of the Hamiltonian the detector is at a; = 0. 



III. DISCRETISATION 

We model the electromagnetic field and the modes of 
the detector with a set of discrete modes. The wavefunc- 
tion of the system can be written as: 

m)) = a\e,0,0) +J2bk\g,lk,0) +J2ck,^\9,0,lk,^) 

k k,uj 

(7) 
where the states involved are product states and, for in- 
stance, \g,lk,0) = |3)|lfe)|0) where \lk) means one pho- 
ton emitted in the fc-th mode and |0) is the zero-quanta 
state of the detector. 

Thus, the initial state vector of the system is |e, 0, 0) 
and the amplitudes obey the Schrodinger equation: 



ia = Oa + ^e^'^^^efcfefc (8) 

k 

(10) 



ICk 



k' ,10 

^Ck,uj + ^ \/rjk^C{k, k')bk 

k' 



where k is the index of the discrete modes used to model 
the electromagnetic field and fc, cu the indexes for the dis- 
crete modes for the k and u of the detector quanta. In 
case k appears by itself, it simply is the value of k of the 
mode. 

Consider for the moment the limit of our Hamiltonian 
that corresponds to the Hamiltonian employed in [3|, i.e 
C{k,k') = Sk,k' and xd — 0. Then the differential equa- 
tions for the amplitudes would be: 



la 



na + ^^kh 

k 

ibk = ^kbk + Cfc" + X! VVkCk.^ 

iCk,LO = ^Ck,i^ + y/Vkbk- 



(11) 

(12) 
(13) 



In this set of equations Ckuj is coupled only to one mode 
of the electromagnetic field, which means that the detec- 
tor interacts immediately with the emitted photon, with- 
out allowing any time delay associated with the distance 
it has to travel from the TLA to the detector. On the 
contrary, in Ea. (|10|l the detector modes interact with a 
superposition which allows for spatial localization of the 
interaction, accounting thus correctly for the time delay 
and the detector position. 

We proceed with a numerical solution for the system of 
differential equations. The discretisation scheme j^, [i3 
is as simple as possible. We choose a range for k and uj 
which we span with equidistant modes. The results are 
considered converged if unaltered upon increasing both 
the range and the density of discrete modes. Of course 
the choice of discretisation range is not arbitrary but 
based on the particulars of the problem. In this case, 
we take |k| around the transition frequency of the TLA 
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FIG. 2: The decay rate of the TLA excited state population, 
coupled to the detector, over the free decay rate, following KS 
approach. Parameters are 27rA/7 = 100,?)/7 = 1, 10 for the 
solid and dashed curve respectively. Together we show, for 
the sake of comparison, the decay rate in the absence of the 
detector (dot-dashed curve). Since we take into account only 
a finite bandwidth, the decay rate is at f = 0. Discretisation 
range: to 20, with 100 modes for the electromagnetic field 
and lOO(k) x 100(ci;) modes for the detector in the same range. 



and the same for lo. Due to the finite interval of |k| 
space that we take into account, the "jump time" is not 
infinitesimal, although it can be made as small as com- 
putationally feasible. In any case, a non-zero t should 
make QZE easier to observe. 

The situation we have considered is equivalent to the 
TLA placed at the center of a (hollow) spherical detector, 
which effectively is a 1-D problem. In this case we have to 
limit to outgoing waves, thus restrict k to positive values. 



IV. RESULTS 

First, we establish a direct correspondence with the 
results obtained in Ref. 1]. We set Xd — 0, C{k,k') — 
6k, k' , the atom-photon coupling independent of k and 
assume that the coupling between the photon and the 
detector is: 



Vk = 



r;/27r 



l + [{k-n)/A\ 



(14) 



with A a measure of the photon energy range for which 
the detector is sensitive and n a parameter defining the 
sharpness of the detector response (n = 6). In Fig. |3we 
show our results which match those obtained in Ref.[l| 
analytically, except for a factor of 10 in the value of rj, 
which we attribute to a possible misprint in the caption 
of their figure; especially since we are unable to repro- 
duce their graph by employing their formula. The initial 
fast drop of the decay rate is due to the finite range of 
frequencies we consider in the discretization, the width 
of this region is of the order of 1/Ai7, where Ail is the 
bandwidth of the discretization. After this transient re- 
gion, the rise of the decay rate to its assymptotic value 
is resolved in accordance Figure 3 of Ref. Q . 

We proceed by considering a detector with finite width, 
the same in all other aspects with the detector in Ref. jj] ■ 
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FIG. 3: The decay rate of the TLA excited state population 
over the free decay rate as its relative position with the detec- 
tor is varied (inset). We considered four cases, from the TLA 
being at the center of the detector (a) to the TLA outside of 
the detector (d). The decay rate changes smoothly from the 
results obtained in Ref.[y| to the free decay rate. The detec- 
tor is assumed to have an effective gaussian profile, with a full 
width half maximum of 33 {h = c — fi = 1), as shown in the 
inset. Parameters: 27rA/7 — 100, rj/'j — 10, discretisation as 
in FigH C{k,k') = 0.103e-«'=-'=')/'5-5)'. 



that the decay process is decelerated, with decay rate 
similar to the one obtained in Ref.0] (^^0.40 vs. ~0.35, 
case (a)). Once the TLA starts to get separated from 
the detector its decay rate approaches fast the free decay 
rate (cases (c) and (d)). The influence of the relative po- 
sition of the detector and the TLA on the dynamics of 
the system is shown in Fig. ^ The time evolution of the 
intensity profile of the emitted photon shows two quali- 
tatively distinct features. In the case where the TLA and 
the detector overlap, the detector captures the emitted 
photon instantly and acts as a "memory" retaining the 
photon close to the TLA and slowing down the decay 
rate. In fact, this is the effect reported in ReL|l|. Once 
the TLA is separated from the detector, the photon trav- 
els uninterrupted until it is absorbed by the detector. In 
this case the atom decays with the free space rate without 
being influenced by the detector. 



V. CONCLUDING REMARKS 
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FIG. 4: The intensity profile of the emitted photon for the 
cases (a) and (d) of Fig|3| In the top part of each figure we 
show the corresponding detector profile and position. The 
distance is measured from the TLA and in units of detector 
full width half maximum. Lighter shading stands for higher 
inensities. 



In Fig. O we show the decay rate of the population of 
the |e) state over the free decay rate, as a function of i in 
cases where the TLA overlaps the detector and where it 
is spatially separated. When they overlap, it is evident 



Starting with the formalism of Ref.[l|, we modified the 
Hamiltonian representing the TLA, electromagnetic field 
and interaction with the detector, so as to explicitly take 
into account the position and the spatial width of the 
detector. The modified formulation allows the analysis 
of the realistic situation in which the detector is spa- 
tially separated from the atom, yielding the model of 
Ref. (l| as a special case which is shown to correspond to 
the detector overlapping with the atom. This is actually 
the case of an excited atom decaying inside a dielectric 
[llL Il3 | . Having calculated the decay probability of the 
TLA, we find that it is not affected by the measurement 
precedure, except in the rather unphysical situation in 
which the atom overlaps with the detector. We are thus 
compelled to conclude that the QZE does not occur _b 
indirect measurements, at least in the context of Ref. 
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